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Abstract 

In a recent paper [M. Colangeli et ai, J. Stat. Mech. P04021, (2011)] it was argued that the Fluctuation 
Relation for the phase space contraction rate A could suitably be extended to non-reversible dissipative 
systems. We strengthen here those arguments, providing analytical and numerical evidence based on the 
properties of a simple irreversible nonequilibrium baker model. We also consider the problem of response, 
showing that the transport coefficients are not affected by the irreversibility of the microscopic dynamics. 
In addition, we prove that a form of detailed balance, hence of equilibrium, holds in the space of relevant 
variables, despite the irreversibility of the phase space dynamics. This corroborates the idea that the same 
stochastic description, which arises from a projection onto a subspace of relevant coordinates, is compatible 
with quite different underlying deterministic dynamics. In other words, the details of the microscopic 
dynamics are largely irrelevant, for what concerns properties such as those concerning the Fluctuation 
Relations, the equilibrium behaviour and the response to perturbations. 
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I. INTRODUCTION 



Physical laws describing the evolution of the microscopic constituents of macroscopic objects 
are commonly assumed to be time reversal invariant, and as such they are usually considered in 
statistical mechanics. We call T-symmetric an evolution which is time reversal invariant. Differ- 
ently, the laws of Thermodynamics, which describe the macroscopic realm, are irreversible. The 
theoretical investigation of the so-called Fluctuation Relations (FRs), which began in the early 
1990's, produced one quantitative approach to emergence of irreversible behavior from the statis- 

Jin, 

tics of dissipative, but time reversal invariant, dynamics [|1L|20. The FRs have then been extended 
and applied to numerous nonequilibrium phenomena ISl-lg]. The usual derivation of the FRs for 
deterministic dynamics relies on the time reversal invariance of the dynamical equations; likewise, 
the derivation of FRs for stochastic dynamics rests on a form of "reversibility" which is relevant 
to the mesoscopic level of description. Stochastic dynamics, thought to represent a projection 
of the phase space dynamics on the space of physical observables (e.g. the 1 -particle space of a 
many particle system) is intrinsically irreversible but, under certain conditions, it is characterized 
by fluctuations which allow every path in the state space, which visits certain states in a given 
chronological order, to be coupled with with another path, which visits the same states in reverse 
order. This property can the be used in the derivation of the FRs for stochastic dynamics. More 
precisely, the T-symmetry present in the stochastic approach to the FRs takes the form of a local 
detailed balance condition, for nonequilibrium systems driven to a steady state, [Q, sj]. 

Recently, the possibility of relaxing the T-symmetry, which is normally assumed to concern 
the whole phase space or the whole state space, have been investigated. In Ref.[9] time reversal 
has been considered in cases in which the dynamics evolves on a lower dimensional manifold, 



llll 



which is not time reversal invariant, but is embedded in the phase space which is. In Refs. [[lo[ 
a similar investigation is performed for stochastic evolutions, in order to identify the minimal 
ingredients required to obtain the FRs or relations of irreversible thermodynamics, such as the 
Onsager-Casimir Relations, previously based on time reversal invariance. 

In a previous paper, [|l2|l . we too considered the question of how fundamental the T-symmetry 
is, in general. To that purpose, we investigated the effect of a "tunable" source of irreversibility on 
a variation of the deterministic and T-reversible toy model known as dissipative Baker Map [|l3|l . 
In this paper, extending the work of [1120, we test the robustness of the FR and of the linear regime, 
with respect to perturbations of the reversibility of the "microscopic" dynamics. This uses some 
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of the insights of the previous work lll2ll . which showed that smoothness of the invariant measure 
and the continuity of the time-reversal symmetry are not necessary in order to derive a FR. The 
conjecture which our data are meant to support, from the point of view of a very idealized setting, 
is that the time evolution of macroscopic observables of real physical systems is compatible with 
many possible underlying dynamics, including irreversible ones. This depends of course on the 
kinds of observables at hand. Hence, we adopted a very simple model in which, however, one may 
distinguish relevant from irrelevant observables, as far as the phase space contraction, denoted by 
A, is concerned. Then, we introduced a source of irreversibility in the dynamics, which affects 
only the irrelevant variable but which affects the structure of the steady state probability distribu- 
tions, including the "equilibrium" ones, and we asked what consequences that may have on the 
validity of the A-FR, namely the FR for the quantity A. 

In our model, the T-symmetry is relaxed by considering a dynamics which is non-invertible. Map- 
pings which are not homeomorphisms constitute a subclass of non-reversible mappings. We expect 
our results to hold, more generally, also in some non-reversible but invertible dynamics. To es- 
tablish a link between deterministic and stochastic dynamical systems, it proves useful to treat 
systems which are non-reversible in the usual deterministic sense (to be specified below), but 
which may be endowed with a weaker notion of reversibility (also to be specified below). For 
example, one may consider a suitable weak notion of reversibility in a projected dynamics (cf. 
Sec. Uni). In particular, our low dimensional toy model may be useful to illustrate how such a 
weakly reversible projected dynamics is consistent with many, possibly irreversible, deterministic 
dynamics. 

Our paper is organized as follows. 

In Secim we describe the general setting and focus on the equilibrium version of the model. The 
proposed distinction between "relevant" and "irrelevant" degrees of freedom, with respect to the 
specific observable A, enable us to verify that the equilibrium version of the A-FR is not affected 
by our source of irreversibility. In SecHlll we derive a detailed balance condition in the projected 
"relevant" space, the one that would correspond to the /i-space of kinetic theory, in the case of 
a system of particles. This notion of detailed balance stems from the more general condition of 
equilibrium in the full phase space, and is proven to work even in presence of irreversible dynam- 
ics. In SeclTVl we consider a nonequilibrium version of our irreversible deterministic map, and 
study the validity of the A-FR, of linear response and of Green-Kubo-like relations for such a map, 
comparing ours with other approaches. Conclusions are drawn in Sec. |Vl 
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Our results can be summarized as follows: 



The condition of detailed balance, in the projected relevant space (e.g. the /x-space in kinetic 
theory), may be derived even from an irreversible equilibrium deterministic dynamics in 
phase space. This shows that a stochastic process derived from a suitable projection onto a 
proper subspace is compatible with many different underlying deterministic dynamics. 

Analytical and numerical results prove the validity of the A-FR for a deterministic dynam- 
ical system which is not time reversal invariant and, yet, satisfies a milder, stochastic-like, 
notion of reversibility. This form of reversibility merely requires the existence of pairs of 
conjugated paths in phase space giving rise to opposite phase space contractions. 

Linear response theory does not hold in our simple models, which violate the conditions re- 
quired by the methods discussed in Refs.||2, 14] and, in particular, do not enjoy any property 



of the kind of Local Thermodynamic Equilibrium. However, the transport coefficients for 
the corresponding reversible and irreversible dynamics coincide. 
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FIG. 1 : The map M in Eq. ([T]) for general values of the parameters £ and q. 
II. GLOBAL CONSERVATIVITY VS. LOCAL DISSIPATIVITY 



Let us introduce the dynamical system {b{,M,n), with phase space U := := M^/Z^ and 
mapping M : U ^ U defined by: 



/ Xn+1 ^ 
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(1) 



for q e [0, and i G (0, ;|], cf. Fig. [T} and with natural measure /i. 



5 



The Jacobian determinant of this map takes the values: 

1 q 



Jm(x) 



Ja 



Jb = 1 



Ai 2i 



for < a; < 



l-2i 



Jc = l + '2q 



JD = 4:i + 2q 



fori<x<- 
2 



tor - < x < - 
2 ~ 4 



for 4 < x < 1 



(2) 



in the four different regions of U. This model generalizes the one introduced in il2|] . as it features 
the two parameters £ and q, which can be tuned to produce different forms of "equilibrium", i.e. 
of natural measures p, which are called non-dissipative steady states because are characterized by 
vanishing phase space contraction rates. In our case, it proves convenient to determine the projec- 
tion of the invariant probability density on the x-coordinate of this map. This can be accomplished 
by integrating over the y-direction the Perron-Frobenius equation [15] for the measures defined on 
the square. This, indeed, yields the evolution equation for the probability measures defined on the 
X axis which are evolved by the map of the interval [0,1] obtained by projecting M on the x axis. 

The calculation of the marginal invariant probability measure can the be performed by intro- 
ducing a Markov partition of the unit square, consisting of two regions: [0, 1/2) and [1/2, 1], 
respectively furnished with the invariant densities pi{x) and Pr{x). Then, the transfer operator 
T associated with the projected dynamics, via the projected Perron-Frobenius equation, can be 
written as: 



/ Pl{Xn+l) 1 




' Pl{Xn) \ 


I PriXn+l) j 




I pr{Xn) ] 



(3) 



where T is defined by: 



T 



(4) 



1 - 2i 1/2 
21 1/2 

The matrix T satisfies the Perron-Frobenius Theorem, hence its largest eigenvalue, A = 1, is sep- 
arated from a spectral gap from its other eigenvalue. Then, the calculation proceeds by evaluating 
the eigenvectors of the transfer operator corresponding to the dominant eigenvalue. The result 
of this procedure shows that the invariant probability density of the map ([T]), projected onto the 
X-axis, depends on the value of £, but not on q, because q only affects the dynamics along the ver- 
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FIG. 2: Result of a numerical simulation for the invariant density of the equilibrium reversible map derived 
from Eq. ([T]) with q = and £ = 0.15, obtained by evolving a set of 2 • 10^ initial conditions randomly (and 
uniformly) chosen on the square [0, 1] x [0, 1]. This shows an invariant density p{x,y) which is uniform 
along the y-coordinate and piecewise constant along the j;-coordinate, attaining the values pi{x) and pr{x) 
for, respectively x G [0, 1/2) and x G [1/2, 1], given in Eq.©. 

tical direction. The corresponding projected density, cf. Fig|2l is given by the piecewise constant 
function: 

2 1 
^'^"^^"TTi^ forO<a;<- 

P{x) = I , (5) 

8i 1 
Prix) = for - < a; < 1 

The marginal probability density p suffices to compute the statistical properties of phase func- 
tions such as the phase space contraction rate A(x, y) = — log J{x, y), because of the special form 
of the Jabobian determinants (|2l), which depend on the x-coordinate only. In this case, one has: 

(A) = — \og J {x,y)fi{dx X dy) = — \og J{x)p{x)dx 
Ju Jo 

The quantity (A) is represented in Figl3]as a function of the parameters {£, q). The figure shows 

that "equilibrium", i.e. by definition the condition in which (A) vanishes, holds only for g = 0, 



independently of the value of C. The dynamical system of Ref. [|12|l can be seen as a special case 
of our map, in which q = \ — 2l and, correspondingly, the only possible equilibrium state of that 
map is given by the further choice ^ = \- 

In this Section, we focus on the case g = and begin by considering the average phase space 
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FIG. 3: The average phase space contraction rate (A), as a function of the parameters £ and q. 




{ 1/2 3/4 1 1/2 1 



FIG. 4: Involution G defined in Eq. Blue lines: the two diagonals, along which the map G reflects the 
two halves of the phase space. 

contraction rate which, in this case, vanishes V£ G (0, |] and can be written as: 

{A) = -\nJAPii-\nJnj = (6) 

The choice i = \, in particular, ensures that the mapping is locally conservative, i.e. that 
A{x,y) = uniformly on U, as all Jacobians are unitary. In particular, i = j leads to a uniform 
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invariant measure /i, which we call "microcanonical" by analogy with the statistical mechanics of 
an isolated particle system with given total energy. The projection of fi is then uniform along the 
X-axis. 

For ^ 7^ ^ and g = 0, the invariant density along the s-axis remains smooth, except at one 
point of discontinuity, x = 1/2 and, in spite of the fluctuations of the phase space volumes, the 
invariant measure is still uniform along the stable manifold, i.e. the vertical direction, as illustrated 
by a numerical simulation reported in FiglU Thus, for i ^ j, we obtain a form of fluctuating 
equilibrium, which we call "canonical" by analogy with the statistical mechanics of a particle 
system in equilibrium with a thermostat at a given temperature ll^ . 

Let us also observe that our equilibrium dynamics (q = 0) are time reversal invariant, according 
to the standard dynamical systems notion of reversibility [|l9ll . because there exists an involution 
G -.U ^lA such that 

MGM = G (7) 

which attains the form: 




G 




f 2x\ 
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2x 



1 \ 



for < a; < I 



for i < X < 1 



(8) 



The mapping G in Eq. ([8]) reflects the half squares [0, 1/2) and [1/2, 1] along the respective 
diagonals, drawn from their lower left to their upper right corners, cf. Fig. |4l 
Hence, according to the definition Eq.©, the baker model in ([T]) with g = is T-symmetric for all 
values of the parameter ^, as shown in Fig. [51 

Consider, now, a trajectory of n time steps, {xq,x^, ...,x„}, along which the average phase 
space contraction rate is given by: 



^ n— 1 

A„(xo) = V In Jm(M' 

n ^ — ^ 



(9) 



fc=0 
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FIG. 5: Check of reversibility for the map ([T]). 
The average phase space contraction rate over the time reversed path is given by: 

n— 1 

A„(GM"(xo)) = V In JMfM'^GMX) 

n ^-^ 

k=0 
^ n—l 

= y In Jm{GMM'^''-^^-%) 

k=0 

Then, the relation Jm{x) = J^liGMx), yields the known result [120: 

A„(G'M%) = -A„(aio) • 

Since the Jacobians (O depend only on the x-coordinate and are piecewise constant, the expres- 
sions © and (flOl ) take the simple form: 

^ n— 1 

A„(io) = y^ynJuW 



n ^-^ 

k=0 
^ n—l 

KXQin-i) = y^lnJM{Qiin-i 



-i)-k 



-A„(zo) 



(11) 
(12) 



k=0 



with ifc the region containing the point M'^Xq, out of the four regions {A, B,C, D} and Q = GM, 
where 

QA = D, QD = A, QB = B, QG = G . (13) 
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Thus, the computation of A„ for the forward (respectively, time reversed) path, can be conveniently 
performed by keeping track only of the coarse-grained sequences of visited regions: 

{ik} = iio,ii,...,in-i) (14) 
{Qi(n-l)-k} = {Qin-i,Qin-2,--,Qio) , (15) 

rather than relying on the more detailed knowledge of the sequence of points {M'^x} and 
{G'M"~'^x} in the phase space. OTII A considerable amount of information, regarding the micro- 



scopic trajectory in the phase space, is lost by passing from the phase space deterministic dynamics 
to the effectively stochastic process arising from the projection of the dynamics onto the x-axis. 
Such a process is described by a Markov jump process, which yields sequences such as those 
of Eqs. (fT4l) and (fT5l) . Nevertheless, this loss of information is irrelevant to compute the phase 
space contraction rate of sets of phase space trajectories. In particular, we may disregard varia- 
tions of the dynamics internal to the single regions, as long as the resulting internal, or "hidden", 
dynamics preserve phase space volumes and do not affect Eqs.(fTTI) and (fT2l) . These observations 
are relevant for the stochastic descriptions of physical phenomena, which are thought to be based 
on reduced (projected) dynamics of phase space deterministic dynamics. Indeed, the projected 
Perron-Frobenius Equation ([3]) is a time-integrated Master equation for the probability densities 
pi and pr over the coarser, projected, state space defined by the chosen Markov partition. These 
considerations are reminiscent of the fact that thermodynamics, to a large extent, does not depend 
on the details of the microscopic dynamics, hence is consistent with many different phase space 
evolutions. This is consequence of the fact that thermodynamics describes the object of interest 
by a few observable quantities, i.e. in a space of reduced dimensionality, which can be seen as 
a projection of the whole phase space. Thus, in order to compute the quantities of interest, one 
may conveniently choose the detailed microscopic dynamics which most easily represent the phe- 
nomenon under investigation. For instance, in our idealized setting, the equilibrium dynamics and 
the existence of some symmetries relating forward and reversed paths in a subspace of the phase 
space, may be investigated by means of the map M with q = 0, which is reversible via the involu- 
tion G depicted in FiglH and by means of the projections of M on the relevant directions. The class 
of dynamics which are equivalent from a given, restricted or projected, standpoint could include 
maps which are not even time reversal invariant. Indeed, one may consider microscopic dynamics 
obtained from map ([T]) introducing an irreversible transformation N which does not contract nor 
expand phase space volumes. This can be simply done by e.g. letting flip the y-coordinates of 
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FIG. 6: The map defined in Eq.(fT6l). wliicli spoils the reversibility of the model 



the phase space points of a vertical strip of width e in the region B: 



Vn+l 




1 - Vv 




for X E [x,x + e] and y E [0, ^] 



for X E [x,x + e] and y E [|, 1] 



(16) 



cf. Figl6]for a graphical representation. 

As pointed out in lll2ll . the composed map K = NM is irreversible because it does not admit 
an inverse. Moreover, the irreversible mechanism of the dynamics described in Eq. (fT6l) gives rise 
to an invariant measure which is fractal along the vertical direction, cf. FigjTl and which is strongly 
at variance with its reversible equilibrium counterpart shown in Figl2l Nevertheless, it is clearly 
seen that Eqs.lfTT]) and (fT2l) still hold true for the map K, despite the irreversible feature of the 
equations of motion. In fact, although the irreversible dynamics K no longer admits an involution 
(hence, Eqs.® and (flOl) can no longer be fulfilled), Eqs. (fTTI) and (fT2l) remain unaltered, because 

maps a point x E i into a point x' E i, and it neither contracts nor expand phase space areas. 
Thus, if we replace the dynamics M with the dynamics K and, accordingly, we take Q = GK, 
the sequences (fT4l) and (fTSi) remain unaltered, since they are invariant under the action of an 
irreversible perturbation of the y-coordinate. 

Now, let uj{io,n — 1) C io and uj{Qin-i,n — 1) C Qin-i denote, respectively, the sets of 
points corresponding to the forward (fT4l) and to the time reversed (fT?)) sequences. Trivially, the 
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FIG. 7: Result of a numerical simulation for the invariant density of the equilibrium irreversible map derived 
from Eqs. ^ and ([T6l ) by setting q = 0, £ = 0.15, x = i and e = ^—£, obtained by evolving a set of 2 • 10^ 
initial conditions randomly (and uniformly) chosen on the square [0, 1] x [0, 1]. The density p{x, y) is not 
smooth along the y-coordinate, which is a signature of the strongly irreversible dynamics given by the map 
N in Eq. (O. 



sets of points corresponding to these symbolic sequences have invariant measures ii{uj{io, n)) and 
li{uj{Qin-i,n)), as in the reversible case. Therefore, in spite of the irreversible modification A^, 
we may say that the dynamics enjoy a form of reversibility which is weaker than the standard 
reversibility in phase space, but which cannot be distinguished from that if observed from the 
stochastic (reduced) viewpoint of the projections on the horizontal direction. 

As a matter of fact, time reversibility is contemplated in stochastic dynamics and amounts 



to the requirement that a sequence of events have positive probability if its reverse does [[170 • 
Because our projected dynamics are not affected by the action of A^, on the level of the stochastic 
description, the phase space reversible dynamics of M are equally stochastically reversible as 
the phase space irreversible dynamics of K. Furthermore, the fact that reshuffles phase space 
points within vertical strips of the square implies that the redistributions of mass due to the phase 
space contraction rate of M and to the rearrangement of phase space volumes produced by N are 
indistinguishable on the projected horizontal space and can be quantified by the same observable 
A. In particular, Eq.© holds for both M and K, hence the statistics of all projected observables 
is the same. Only in the case that one is interested on observables which explicitly concern the y 
coordinate would the two dynamics be distinguishable, but as long as one focuses on quantities 
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which do not depend on y, or which result from a projection on the x axis, M and K lead to the 
same conclusions: the corresponding reduced stochastic evolution is exactly the same, as in the 
case of felt dynamics introduced in [|l8l] . 

Therefore, we may state that the map K with g = represents a kind of equilibrium but 
irreversible dynamics. Although this might appear contradictory, it is simply explained by the 
observation that the equilibrium behavior concerns the level of the horizontal projection, which is 
stochastically reversible, while the irreversibility concerns the phase space. This situation differs 
from that of in which time reversible nonequilibrium systems are recognized to be com- 
mon -see, e.g. the standard models of nonequilibrium molecular dynamics- while equilibrium 
irreversible systems are thought to be rare, as far as phase space is concerned. 

Our study concerns, instead, the bridge between deterministic and stochastic-like descriptions. 
In particular, we are going to show that the s-projection of K, being an equilibrium model, satisfies 
the principle of detailed balance (DB), from the point of view of the stochastic dynamics, in spite 



of its irreversibility in the phase space, p 81] 



This naturally connects with the distinction between relevant and irrelevant coordinates which 
underlies the statistical mechanics reduction of deterministic descriptions in the phase space to 
stochastic descriptions, which typically concern the one-particle space 11221 12311. Clearly, our dy- 
namical system is too simple to allow a physically meaningful distinction between relevant and 
irrelevant variables. Therefore, we merely focus on one of them (the x-coordinate) and regard 
the other as irrelevant (the y-coordinate), without implying that one subspace of our subspace is 
endowed with any special meaning. 



III. FROM PHASE SPACE TO DETAILED BALANCE 



In this Section we introduce the notion of detailed balance in the phase space (PSDB) which 
constitutes a strong concept of equilibrium dynamics. We will show that the standard notion of DB 
descends from PSDB through a projection of the phase space dynamics onto a suitable subspace. 
It will also be shown that the DB condition is insensitive to the reversibility properties of the full 
phase space equilibrium dynamics. Let 

x„+i = Mx^ , xeU 
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be the microscopic dynamics, where M is time reversal invariant, with involution G. For sake of 
simplicity, we deal with discrete time, t E Z, but flows S**, t G M could be treated similarly. 

Consider two sets in phase space, U,V CU. Let W = M^U f}V = {xeV : M'^s G U] 
be the set of final points of trajectory segments starting in U, which fall in V after k iterations of 
M, and let M^'^W be the corresponding set of initial conditions. Take an invariant measure fi for 
M, so that niW) = iJ,{M~''W), and let Afc(x) be the phase space contraction along the trajectory 
starting in x G M~^W . Its time reverse trajectory, which gives rise to the opposite phase space 
contraction, cf. Eqs.([TOl) and starts in x = GMM^x G GMW . CaU 

W = M^GMW = M'^GM {M^U nV) = GMU n M^GMV 

the set of final conditions of all trajectories which are time reverses of those ending in W , cf. Fig. 
[HI where the second equality comes from the definition of W and the third equality follows from 
the property M^G = GM~'' of the involution G. The measure of this set is given by: 

H{W) = niM'^GMV nGMU) 

= n{M''M-^GV n M-^GU) = i2{M^-\GV n M-''GU)) 
= ii{M^-\GV n GM^U)) = fi{GV n GM'^U) 

= i^{G{V n M'^U)) = fi{GW) (17) 



We define the phase space detailed balance (PSDB) as the condition for which the probability 
of having opposite phase space contractions are equal: 

^i{M-^W) = fi{M-^W) i.e. n{W) = n{W) (18) 

Becasue of Eq.(fT7l). this condition may also be written as 

jj{W) = i2{GW) (19) 

which is to say that PSDB requires the M-invariant measure /i to be also G-invariant, since W may 
be any subset of W. Calling a function $ odd with respect to time time reversal if = — 

for all X E U, and defining equilibrium the situation in which the mean value of all such odd 
observables vanishes, we obtain that Eq.(fT9l) implies equilibrium. Indeed, take any $ which is odd 
with respect to the time reversal and choose W = {x E U : ^{x) > 0}, so that GW = {x E U : 
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FIG. 8: Set of points, in phase phase, belonging to the sets U and V in the forward (upper path) and in the 
time reversed (lower path) trajectories 



where Jg = 1. In particular, the PSDB implies that the average of the phase space contraction rate 
A{x) = log Jpix)^^ vanishes, as required for equilibrium in SecHlIl 

Is there any relation between PSDB and the standard DB, which implies equilibrium on the 
level of the projected dynamics (the caricature of the one particle or yU-space)?[|39!] To derive stan- 
dard DB, let us eliminate the irrelevant coordinates, which only contribute to noise, by projecting 
Eq.(fT9l) on the subspace of relevant coordinates, which we call yU-space. This can be done for sets 
of the form W = W^^^ x C„_i and W = W^^^ x Cn-i, where W^^'^ and W^^^ denote the pro- 
jections of the sets W and W onto the yU-space and C„_i, C„_i span the remaining, noisy, space, 
cf. the hyper-cylinders illustrated in Fig. |9l where W^^'^ and W^^'^ denote, respectively, the event 
on the /i-space. Then, if we denote by fi^^^ the measure in the /i-space induced from the invariant 



$(x) < 0} and = for all x e W \ U GW). Then the following holds: 




(20) 
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FIG. 9: The sets VF^^-* and GW^^\ after reduction from phase space to the ^-space 



measure on the phase space, we have that: 

iw^^)) = I ^{dx) , = / fi{dx) (21) 

Detailed balance holds if 

^(1) (1^(1)) = ^(1) , or (1^(1)) = {GW^'^) (22) 

Thus, PSDB implies the standard DB, since DB amounts to the condition of PSBD restricted to 
sets W and W of the form here introduced. Moreover, the projection procedure typically smoothes 
out singularities, hence the induced invariant measure /x*^^^ usually is regular and has an invariant 
density p^^^.H 

In the derivation of the stochastic description as a projection of some deterministic phase space 
dynamics, the DB condition (|22l) is usually assumed to be the consequence of the time reversal 
invariance of the microscopic dynamics, once equilibrium is reached. In our investigation this 
amounts to require the existence of the involution G, defined in the phase space. We are now 
going to see that this requirement may be relaxed in simple cases, such as those of the irreversible 
dynamics K discussed above. Using the same notation, Eq.(fT9l) can then be written as: 

/i(Mi n j) = iJ,{MGMj n GMi) (23) 

where jj,{Mi fl j) is the conditional probability of being in the region j one time step after having 
been in region i, with i, j E {A, B, C, D}, our finite state space. The quantity /x(Mz fl j) may be 
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p 



1 1\ 

2 2 



(24) 



rewritten as iJi{Mi fl j) = p{j\i) = Pij, with notation reminiscent of stochastic descriptions. The 
Pij's then constitute the elements of the transition matrix: 

/ 

2£ 1 - 2£ 
H 
\ 2£ 1 - 2£ / 

which defines a j'toc/ia^ric process for the dynamics of the state in {A, B,C, D}. Application of the 
Perron-Frobenius Theorem to the matrix P reveals the existence of one left eigenvector of P, as- 
sociated with the eigenvalue A = 1 (i.e. the geometric multiplicity of such eigenvalue is 1), which, 
thus, implies the existence of a unique (coarse-grained) steady state measure (/i^, /ie, /ic, 
where: 

f 91 

if i = A,C,D; 

(25) 

if i = B. 



2e 



Eq. (l25l) highlights the fact that the measure in the full phase space is uniform along the stable 
manifold and piecewise constant along the unstable one, as previously illustrated in Figl2l This 
comes from the fact that, since the regions {A,B,C,D} are determined only by the x-coordinate, 
the stochastic transition matrix P does not depend on the dynamics along the stable manifold, 
and, hence, = ( J")^^. As a result, the invariant measure in Eq. (l25l) depends only on x and is 
constant on y. 

Consider, for instance, the one-step transition A ^ C, whose probability is the measure of 
W = MA n C, which equals that of the set of the corresponding initial conditions M~^W = 
A n M~^C. The probability of the reverse transition QC = C — t- QA = D, where we have 
recalled the relations (fT3l) . is the measure of the set W = MQC fl QA, cf. Fig{TOl Since the sets 
W and W span the whole range [0, 1] in the vertical direction, as in Fig. |9l the measures fi{W) 
and n{W) can be calculated just in terms of the transition probabilities (|24l) and of the projected 
invariant measures (|25]) . The result is 



fi{MA nC)= fiAPAc 



l^cPcD = KMC n D) 



(26) 



1 + M 

for any L Hence, DB holds as expected, because we are dealing with an equilibrium case, al- 
though the underlying dynamics is irreversible. It is now interesting check what happens when the 
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M QCnQA 



MA n C 



FIG. 10: Sets undergoing the forward path {A — > C, blue arrow) and the time reversed path {QC — > QA 
black arrow) for the map ([T]) with q = 0. 



microscopic dynamics is pulled out of equilibrium. For instance, consider K = NM, with as 



in (fT6l ) and M the time reversible mapping ([T]), with q 



2i. For this map M, we may consider 



an involution G which is consistent with the following equalities (cf. Eqs.(36) in [12]): 



QA = A, QB = C , QC = B , QD = D 



(27) 



hence which differs from the G in Eq.®. Then, the time reverse of the transition A — )■ C is given 
by QC QA i.e. B ^ A, and we get: 



fi{MA nC) = haPac 



l + 4i 
(1 - 2i)^ 



(28) 
(29) 



which shows that (l28l) and (|29l) do not coincide, and that PSDB and DB are violated for q ^ 0, i.e. 
outside the equilibrium defined via the chosen G. In this case, only i = j leads to the equilibrium 
state which is, in addition, microcanonical ll4in . 



IV. THE FLUCTUATION RELATION AND NONEQUILIBRIUM RESPONSE 



The fluctuation relation for A, the A-FR, originally proposed by Evans, Cohen and Morriss 
11241. and developed by Gallavotti and Cohen concerns the statistics of the mean phase space 
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contraction rate A„ , over the steady state ensemble of phase space trajectory segments of a large 
number of steps, n. Equivalently, it concerns the statistics of A„, computed over segments of a 
unique steady state phase space trajectory, broken in segments {xi, x„} of length n. 
The dynamics are called dissipative if 

(A) = / A{x) fi{dx) > 
Ju 

where (A) is the steady state mean of A, i.e. it is computed with respect to the natural measure fx 
on U. It is convenienet to introduce the dimensionless phase space contraction rate e„ = A„/ (A) 
because its range, [— Amax/ (A), Amax/ (A)], does not change with n, while the values taken by 
e„ tend to more and more densely fill it as n grows. Further, we denote by nn{Bp^s) the probability 
that e„, computed over a segment of n-steps of a typical trajectory, falls in the interval Bp^s = 
(p — S,p + S), for some fixed 5 > 0. In other words, one may write 

-KniBp^s) = fJ.{u} A,n) , whcrc wa,„ = {x G W : e„ G Bp^s} (30) 



For growing n, iXn peaks around the mean value (e„) = 1, but fluctuations about this mean 



occur with positive probability at any finite n. In particular, under certain conditions, [lll,|5|. 
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may 



TCn obeys a large deviation principle with a given rate functional (, in the sense that the limit 

lim TTniBp^s) = e-'^K(^')+^^] (31) 

exists, with < S. In particular, if the support of the invariant measure is the whole phase space 
U, time reversibility guarantees that the support of 7r„ is symmetric around 0, and one can consider 
the ratio 

'n'njBp^s) 

T^niB^p^s) 

In our case, this ratio equals the ratio of the measures of a pair of sets conjugated by time reversal, 
as in Eqs. (fTTI) and (fT2l) . Then, the validity of the A-FR means that there exists p* > such that 

r ^ v 1 1 ■■ en{x) e Bp^s}) . . r 
p-6< hm — — log — . . „ ^P + ^ (32) 

if \p\ < p* and 6 > 0. 

So far, the proofs of this and other FR's appeared in the literature, notably those for the fluctua- 



tions of the Dissipation Function fl, ^ 



21\\ . rely on the existence of an involution representing 



time reversal in phase space, while they rely on the principle of microscopic reversibility in the 



20 



state space of stochastic processes. So, whatever the context, the relevant notion of time reversibil- 
ity has always been used. Therefore, if time reversibility is broken, but is broken as in the case of 
the map K, which enjoys a weaker form of reversibility requiring only the existence of the pairs 
of conjugate trajectories (fT4l) and (fT5l) . the A-FR should remain valid. Indeed, the existence of 
this weaker reversibility is consistent with the principle of microscopic reversibility, adopted in 
the stochastic approach by e.g. Lebowitz and Spohn, [81]. 

Let us then investigate the validity of the A-FR for the deterministic model K = NM with 
given by Eq.® and M by Eq.O, which may or may not lead to equilibrium, depending on 
the value of the parameter q. As discussed in Sec. |nl is irreversible, M is reversible and the 
mapping Q = GK, appearing in the definition of the time reversed path (fT5l) . is properly defined 
in both cases. Take q = \ — 2l, with I ^ \, consistently with [12]. Then, the A-FR may be written 
as: 

p-^<lim4^1os ;'"''°-""'", <P + ^ (33) 
with \p\ < p*, for some p* > and any 5 > 0. To prove Eq.(l33]) for the map K = NM, one must 
compute the invariant probability measure 7r„ and let n grow without bounds. 



This is guaranteed by the proof of the validity of the A-FR given for M in Ref.ftl^]. which only 
relies on the invariant measure in the projected space. 

We illustrate this result by means of numerical simulations, which we have performed for 
different values of n, with £ = 0.15, cf. Fig{TTJ The numerical simulations, shown in Fig. [121 
show how the large deviation rate functional C, is generated, and that it is smooth and strictly 
convex in the whole range of observed fluctuations, as required by the theory of the A-FR. In 
Fig{T3] we also plotted the expression in the center of Eq.(l33l). which is consistent with the validity 
of the A-FR. Once the A-FR is proven to hold for the map K, one may be tempted to assess the 
validity of the Green-Kubo formulas as well as of the Onsager reciprocal relations, by following 
e.g. the strategies of HI in the limit of small external drivings, as summarized in [|3|,|4|]. 

To this end, we consider A as the entropy production rate. [|42n and we briefly summarize the the 
argument, for sake of completeness. 



The main steps are the following 



Is, 



assume that the system is subjected to k fields F = (Fi, F2, .., Fk), that A vanishes when all 
drivings vanish and that 

k 

A(x) = 5^F,J°(x) + 0(F2) (34) 



e=i 
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FIG. 11: Probability measure iTn{Bp^s) 
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FIG. 12: Rate functional ^„ associated with vr^for different values of n. As expected from theory, the 
curves Cn are expected to move downwards for growing n, so that, in the n — > oo limit, ({p) intersects the 
horizontal axis only inp = 1. 



which defines the currents J°, which are proportional to the forces F^. 
• The decay of the A autocorrelation function required for the A-FR to hold, leads to the 
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0.8 



FIG. 13: Check of the A-FR, Eq. (1331 ) for the irreversible dynamics dictated by the map K = NM. Blue 
points: results of the numerical simulations with n = 2 • 10^. 



following expansion for the rate function C,: 

Cip) = - 1)' + Oiip - ifF') (35) 

where C2 is related to the time autocorrelation of A. In other words, the rate functional 
is quadratic for small deviations from the mean p = 1, in accord with the Central Limit 
Theorem. 

• Introduce the nonlinear currents as Je{x) = 9ir^A(x), and the transport coefficients = 

{Je)\F=o- Then, one obtains: 

(A) = ^ 5Z idpAJi) + dFMr))\F=oFiFr = ^ ^^'^ + Lre)FeFr (36) 

£,r=l e,r=l 

to second order in the forces. 

• Equation (|33l) implies (A) = Thus, equating the latter expression with Eq. (l36l) and by 
considering (L^^ + Lri)/2 with i = r, one obtains the Green-Kubo relations. 

In our case, the rate functional is clearly quadratic, as shown by our simulations of the dynam- 
ics of K. In particular, the red quadratic curve in Fig [14] reproduces nicely the behavior of the 
numerical data for the rate functional Cn> corresponding to trajectory segments of n = 200 steps. 
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FIG. 14: Rate functional associated to vr„(i3p 5). Blue points: results of the numerical simulation with 
n = 10^. Red line: Fitting of numerical data with the parabola a(a;— 1)^+6, with parameters a = 2.66-10"^ 
and b = 2.68 • 10~2_ 



The necessity for a parameter 6 7^ in the parabola is due the finiteness of n: indeed 6—7-0 when 
— 7- 00. However, in spite of the validity of the A-FR for the irreversible map K, which entails 
that the irreversible map behaves to some extent equivalently to the reversible map M, the argu- 
ment of [|91] leading to the Green-Kubo relations cannot be reproduced here. In fact, it relies on the 
differentiability of the SRB measure as well as on the reversibility of the microscopic dynamics, 
which are both violated in the case of K. Alternatively, one may think of deriving linear response 
from the A-FR through the approach of Ref.[|20 (SRE, hereafter), which does not explicitly require 
the differentiability of the invariant measure and the reversibility of the dynamics. In particular, 
SRE deals with a Nose-Hoover thermostatted A^-particle system and obtains: 



FeV 

2Ko 



2L{F,)kBT 
V 



+ 



nN 



(37) 



for the variance cr|- of A„, provided A can be identified with the dissipation function Vt. Here, 
where Kq is the target kinetic energy for the thermostatted particles, which corresponds to the 
inverse temperature /3, is the number of particles, V is the volume and F^. denotes the external 
force (i.e. the bias) acting on the system. The quantity L{Fe) is defined by 



(38) 



L(Fe) = / dt{{^{t) - (vI/))(vI/(0) - m) 
24 



and 1/(0) = limp^^o L(Fe) is the linear transport coefficient. The derivation of the Green-Kubo 
formulae is completed by comparing (|37l) with the relation (A) = {'^)Ff> = |(t|- (Fe), which is 
implied by the FR, which yields: 



L(0) = lim 



(VI/) 



PV / t/t(vl/(t)vl/(0)) 

Fe^O he Jo 

In our framework of simple dynamical systems, the "current" vp could be defined as: 

forxe A,D 

= ^ 1 forx G B , 



for X E C 



(39) 



(40) 



which implies an average current {'^) = (1 — 4£) /(I + 4£), cf. Eq.(41) in Ref. [|l2h . where the role 
of the external force Fe was played by the bias 6 = 2 — 1/(1 — 2£). ll43ll 

Nevertheless, following SRE may be problematic, as N is required to be large enough in order 
to derive (|39l ) from (ITTI) . something which cannot be granted in low-dimensional systems as the 
map under consideration. Indeed, our numerical simulations reveal that an interesting scenario 
arises in the computation of the quantity L{Fe), which may be conveniently approximated by: 



L{Fe 



Ne. 



(41) 



A:=0 



where the upper limit of the integral in (|38l) is replaced by Nuer and the correlations are computed 



over an ensemble of fully decorrelated initial conditions {x^^}, with j 



Nens, picked at 



random so that they occur in the ensemble with the frequency corresponding to the natural invariant 
measure of the dynamical system. Then, the coefficient L(0) may be computed from Eq. (|4ll) 
by considering the limit of vanishing bias, i.e. by taking (vp) = and a microcanonical ||44 1 
equilibrium ensemble of initial conditions {x^ ^}. The values Nens and Nuer must be chosen with 
care, in order to guarantee the convergence of the sums in (|4TI) . cf. FigfTSl 

Our simulations show that the irreversible character of the dynamics, which affects the "irrel- 
evant" variable y, has no influence on the convergence of the Green-Kubo formula Eq. (|4T)) . cf. 
FigUa 

To study the linear response for the irreversible map, that is the existence of the limit L(0) = 
HmiT^^o L{Fe), we varied the value of bias within four different windows of magnitude, cf. FiglTTl 
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FIG. 15: Numerical computation of L(0) from Eq. (|4TI) . for the irreversible map K = MN derived from 
Eqs. (fT6l) and ([U, with q = ^ — 2£, for different values of N^ns- 
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FIG. 16: Computation of L(0) for the irreversible, K = MN, and the corresponding reversible map, 
K = M, with Nens = 5 • 10^. The irreversibility does not affect the convergence in Eq.(l4TI). 

and obtained that, in spite of the validity of the A-FR, no linear response can be claimed for our 
low-dimensional system, unless this is verified at exceedingly small bias. This fact cannot be 
blamed on the irreversible nature of the evolution, since we have observed that the response of the 
irreversible map coincides with the response of the corresponding T-symmetric one, cf. Fig{T8l It is 
more related to the irregularity typical of transport phenomena in low dimension i34l1 . Therefore, 
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FIG. 17: Numerical computation of the coefficient L{Fe) as a function of the bias Fg at different orders 
of magnitude of the bias. The coefficient L{Fe) presents a very irregular structure, which may prevent the 
existence of the limit L(0) = limiT^^^o L{Fe). 



the dynamics of the map K = MN proves that its irreversible component, the map A^, affects 
neither the validity of the A-FR nor the response of the system to an external bias. 

Per se, the fact that neither the approach of |^ nor that of SRE are applicable does not imply 
that no linear response can be established. However, this a clear observation for our model, which 
does not enjoy many properties of the systems of l^and many others of the systems of 1'^. More 
importantly, as pointed out in e.g. Refs.l 16. 35], the physical linear response relies on the 
occurrence of local equilibrium, in the sense that real space may be thought of as a "collection" of 
cells, each of which contains a statistically significant number of interacting particles. Clearly, our 
two-dimensional dynamical system may mimic only a few features of a real A^-particle system, 
and the local equilibrium property is out of question. 



V. CONCLUSIONS 

In this paper, we have studied a kind of baker model whose properties are determined by two 
parameters, one of which, q, may be conveniently tuned in order to fix the distance from the 



"equilibrium" state. In particular, we have generalized the nonequilibrium map introduced in [|l2|l . 
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FIG. 18: Transport coefficient L{Fe) for the irreversible, K = MN, and the corresponding reversible map, 
K = M, with Nens = 5 • 10^. The irreversible component of the dynamics, A^, leaves the response of the 
system unaffected. 

which can be recovered from our present model by a suitable choice of q. The system studied 
in this work is only a caricature of a real particle system subjected to the action of an external 
driving, but can be studied in detail and thus is useful in understanding the projection procedures 
typically thought be necessary to obtain a coarser, stochastic-like, description from a microscopic, 
deterministic, one. 

If one is interested in one specific observable, e.g. A, our model allows a simple identification 
of the relevant and of the irrelevant variables: the phase function A depends on the Jacobians of 
the mapping, and the Jacobians only depend on one of the variables, x, which is then the only 
relevant variable. This implies that adding a source of irreversibility which affects the "irrelevant" 
degree of freedom y, has no effect on the equilibrium state defined by (A) = 0, and on the validity 
of the A-FR. 

Projecting the invariant measure on the reduced space of the relevant variable not only pro- 
duces a probability density, which is smooth along the unstable manifold, except for one point of 
discontinuity, but also shows that equilibrium survives in the reduced space of the x-variable, even 
in presence of irreversible full phase space dynamics. More precisely, equilibrium survives in the 
form of detailed balance, which is the notion characterizing the equilibrium states in the spaces of 
relevant observables, e.g. the yU-space. This is due to the fact that in these spaces the details of the 
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dynamics of the "noisy" degrees of freedom are irrelevant. 

We have also considered the validity of the A-FR and of the transport properties of an irre- 
versible dynamical system. So far the validity of the A-FR has been derived as a property of 
T-symmetric dynamical systems. Our analysis, also supported by numerical tests, shows that as 
long as the source of irreversibility affects only the "irrelevant" degrees of freedom, theA-FR and 
the transport laws (linear and nonlinear response) hold indistinguishably for both reversible and 
irreversible phase space dynamics. Thus, our results extend to phase space dynamics some of the 
considerations raised e.g. in Ref.[l 1] for stochastic dynamics. This is, in fact, done by relating the 
phase space dynamics to its projections which, according to basic tenets of statistical mechanics, 
should result in stochastic evolutions. 

In our investigation on a prototype of irreversible dynamical system, we discussed a model which 

y one way to accomplish that. However, one 



is, manifestly, non Anosov: non-invertibility is mere ^ 
may think of more realistic non-invertible dynamics 113611 . In that case, the orientations of the phase 
space regions in which the phase space is decomposed will play an important role. Nevertheless, in 
more physical, higher dimensional, models, it might not be so straightforward to identify relevant 
or irrelevant variables with fixed directions in phase space, as was possible to do with our model. 
Investigations of further, more elaborate models is, therefore, worthwhile. 
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[37] The last regions in each sequence, i„ and Qi-i, need not be taken into account, as they are unessential 
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[38] Note that DB is often referred to as the principle of microscopic reversibility 11201] . which is then to be 
understood as a notion of reversibility in the reduced space, not in the full phase space. 

[39] Observe that some authors distinguish DB dynamics from detailed balance steady state, which proves 
to be a convenient tool in the analysis of stochastic processes [7]. In this case, a given evolution law is 
called DB dynamics if its steady state is a DB state. 

[40] For instance, the map ([T) with h = ^ — 2i which, for arbitrary i, is dissipative but still equipped with 



a projected invariant density which is smooth along the unstable manifold, 11121] . 
[41] The map featuring q = ^ — 2£ attains equilibrium for g = 0, which gives (■ = \, corresponding, as 
discussed in Ref.|ll20, to a microcanonical equilibrium distribution. 
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[42] This identification must be done cum grano salis, as explained in e.g. ||2|, 
[43] Adopting the definition of current in Eq. ( |40l ) for the map of Sec. JIIJ one finds a "positive current" for 
expanding phase space volumes. This is different from the case of the most common deterministically 
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thermostatted dynamics, typical of nonequilibrium molecular dynamics. However, there is no general 
principle which imposes phase space contraction in presence of positive currents. Indeed, one may 
consider models without any phase space variations, or models with phase space expansion for positive 
currents (see e.g. [33] for certain parameter values). Moreover, in our highly idealized model, the 
observable ^ can be defined differently, so that it may take whatever values one likes. 
[44] See footnote at the end of Sec Hill 
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